Cosmological consequences of a Chaplygin gas dark energy 
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A combination of recent observational results has given rise to what is currently known as the 
dark energy problem. Although several possible candidates have been extensively discussed in the 
literature to date the nature of this dark energy component is not well understood at present. In 
this paper we investigate some cosmological implications of another dark energy candidate: an 
exotic fluid known as the Chaplygin gas, which is characterized by an equation of state p — —A/p, 
where y4 is a positive constant. By assuming a flat scenario driven by non-relativistic matter plus a 
Chaplygin gas dark energy we study the influence of such a component on the statistical properties 
of gravitational lenses. A comparison between the predicted age of the universe and the latest age 
estimates of globular clusters is also included and the results briefly discussed. In general, we find 
that the behavior of this class of models may be interpreted as an intermediary case between the 
standard and ACDM scenarios. 

PACS numbers: 98.80.Es; 95.35.+d; 98.62.Sb 



I. INTRODUCTION 

From a large number of observational evidence, the 
currently favoured cosmological model is flat, accelerated 
and composed of '--^ 1/3 of matter (barionic + dark) and 
~ 2/3 of a negative-pressure dark component, usually 
named dark energy or "quintessence" . The nature of 
such an unclustered dark energy component, however, is 
not very well understood at present, giving rise to many 
theoretical speculations. 

Certainly, the most extensively studied explanation for 
this dark energy problem is the vaccum energy density 
or cosmological constant (A) although other interesting 
possibilities are also alive in the current literature. Some 
examples are: a very light scalar field (f>, whose effec- 
tive potential V{(l)) leads to an accelerated phase at the 
late stages of the Universe [0, a X-matter component 
1^ , which is simply characterized by an equation of state 
Px = ^xPxi where —1<<jJx<Q and that includes, as 
a particular case, models with a cosmological constant 
(ACDM), a vaccum decaying energy density or a time 
varying A-term whose the present value of the cosmo- 
logical constant (Aq) is a remnant of the primordial in- 
flationary/deflationary stage 1^, geometrical effects from 
extra dimensions [Q or still an exotic fluid, the so-called 
Chaplygin gas, whose equation of state is given by 
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where A is a positive constant . 

All the above mentioned candidates for quintessence 
have interesting features that make them at some level 
compatible with the recent obervational facts (see, for ex- 
ample, ^, 1^, H ^). Although most of these scenarios have 
been extensively explored in the recent literature, in the 
case of a Chaplygin gas-type dark energy, however, only 
few analysis have focused attention on its cosmological 
consequences. From a theoretical viewpoint, an interest- 
ing connection between the Chaplygin gas equation of 
state and String theory has been identified |ri|, p^ . 
As explained in ml 
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a Chaplygin gas-type equation 
of state is associated with the parametrization invariant 
Nambu-Goto d-brane action in a d -|- 2 spacetime. In the 
light-cone parametrization, such an action reduces itself 
to the action of a Newtonian fluid which obeys Eq. (1) so 
that the Chaplygin gas corresponds effectively to a gas 
of d-branes in a d -|- 2 spacetime. Moreover, the Chaply- 
gin gas is the only gas known to admit supersymmetric 
generalization pi] ]. 

From the observational viewpoint, it has been argued 
that the Chaplygin gas may unify the cold dark matter 
and the dark energy scenarios 1 12 . The reason for such a 



belief is the general behaviour of the Chaplygin gas equa- 
tion of state: it can behave as cold dark matter at small 
scales and as a negative-pressure dark energy component 
at large scales. Recently, Fabris et al. |1^ analysed a 
cold dark matter plus a Chaplygin gas scenario in the 
light of type la supernovae data (SNe la). As a general 
result, they found a universe completely dominated by 
the Chaplygin gas as the best fit model. More recently, 
Avelino et al. |Q used a larger sample of SNe la and 
the shape of the matter power spectrum to show that 
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such data restrict the model to a behaviour that closely 
matches that of a ACDM models while Bento et al. ||l7|| 
showed that the location of the CMB peaks imposes tight 
constraints on the free parameters of the model. 

The aim of this paper is to explore some other observa- 
tional consequences of a Chaplygin gas dark energy. We 
mainly focus our attention on the constraints from statis- 
tical properties of gravitationally lensed quasars (QSOs) 
on the Eq. (1). We also investigate other observational 
quantities like the deceleration parameter, the accelera- 
tion redshift and the expanding age of the the universe. 
To obtain such results we assume a flat model driven by 
non-relativistic matter plus a Chaplygin gas dark energy 
component (from now on CgCDM). 

This paper is organized in the following way. In Sec. II 
the field equations and distance formulas are presented. 
We also derive the expression for the deceleration pa- 
rameter and discuss the redshift at which the accelerated 
expansion begins. The predicted age of the Universe in 
the context of CgCDM models is briefly discussed in Sec. 
III. We then proceed to analyse the constraints from lens- 
ing statistics on these scenarios in In Sec. IV. We end the 
paper by summarizing the main results in the conclusion 
section. 



II. FIELD EQUATIONS, DECELERATION 
PARAMETER AND DISTANCE FORMULAS 
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FIG. 1: Deceleration parameter as a function of redshift for 
some selected values of Qm and As. The horizontal line la- 
beled decelerating/accelerating {qo — 0) divides models with 
a decelerating or accelerating expansion at a given redshift. 



Let us now consider the Friedmann- Robertson- Walker 
(FRW) line element (c = 1) 



ds^ = dt"^ - R^{t) 



1 — k?'' 



, (2) 



where fc = 0, ±1 is the curvature parameter of the spatial 
section, r, 9, and (j) are dimensionless comoving coordi- 
nates, and R{t) is the scale factor. Since the two com- 
ponents (nonrelativistic matter and Chaplygin gas) are 
separately conserved, we use the energy conservation law 
together with Eq. (1) to flnd the following expression for 
the Chaplygin gas density 



PCg 




or, equivalently, 



PCg = pCgA As + As 



Ro 



(3) 



(4) 



where the subscript o denotes present day quantities, 
B = p^g^ — A and A^ = A/p^^^ is a quantity related 
with the sound speed for the Chaplygin gas today. As 
can be seen from Eq. (|3|), the Chaplygin gas interpo- 
lates between non-relativistic matter {pcg{R ^ 0) ~ 
\/B / R^) and negative-pressure dark component regimes 
{pcgiR'^'x) ~ VA). 



The Friedmann's equation for the kind of models we 
are considering is 



R 



Q^,i^)^ + nCg^As + il-As)i^}^ 



(5) 



In the above equation, an overdot denotes derivative with 
respect to time, Ho = lOO/i km.s^^Mpc^^ is the present 
day value of the Hubble parameter, and rim and ^Cg £^re, 
respectively, the matter and the Chaplygin gas density 
parameters. 

The deceleration parameter, usually deflned as qo = 
—RR/ R^\t„, now takes the following form 



l[nn, + ncg{i~As)] 



Cg 



(6) 



As one may check, for Ag — and Ag — I, the above 
expressions reduce to the standard and ACDM models, 
respectively. 

Figure 1 shows the behavior of the deceleration param- 
eter as a function of redshift for selected values of flm and 
Ag. The best flt ACDM case is also showed for the sake 
of comparison {Ag = 1). Note that the value of Ag de- 
termines the acceleration redshift Za- At late times, a 
CgCDM model with = 0.1 and Ag = 0.9 accelerates 
faster than a ACDM scenario with = 0.3. In such 
a model the accelerated expansion begins at Za — 0.51. 
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FIG. 2: Dimensionless angular diameter distance as a function 
of the source redshift [zs) for some selected values of As. In 
all curves the value of the matter density parameter has been 
fixed (Oin = 0.3). 



For the best fit model found in Ref. [g^, i.e., l^m = 
and As = 0.92, the universe is strongly accelerated to- 
day with the accelerated phase begining at Za — 0.68 
whereas for = 1 and = 0.6 and fim = 0.3 we find, 
respectively, Za — 0.67 and Za — 0.07. 

From Eqs. (2) and (5), it is straightforward to show 
that the comoving distance ri{z) to a light source located 
at r = ri and t = ti and observed at r = and t = to 
can be written as 



(1 + zs) 



(9) 



dx 



x'g f ip^T^an ^s) 



In Fig. 2 we show the dimensionless angular diam- 
eter distance between an observer O and the source S 
[DosHo) as a function of the source redshift (zg) for 
il,n — 0.3 and selected values of As- As physically ex- 
pected, the larger the value of Ag the larger the distance 
that is predicted between two redshifts. This result shows 
that, for the value of f^m considered, the behaviour of 
this class of CgCDM models may be interpreted as an 
intermediary case between the ACDM {As = 1) and the 
Einstein-de Sitter {As = 0) scenarios. This particular 
feature of CgCDM models may be important for the lens- 
ing statistics analysis because, as is well known, the large 
distances predicted by ACDM models make the lensing 
constraints on the vaccum energy contribution very re- 
strictive (see, for instance, ||l^). In this concern, we ex- 
pect that the constraints from this particular test will 
be weaker for CgCDM scenarios than for their ACDM 
counterparts. It is worth mentioning that the behav- 
ior of CgCDM cosmology can be very different from that 
one present by ACDM scenarios and general quintessence 
cosmologies. For example, as shown in Ref. [|l9[, the tra- 
jectories of the statefinder parameters j20j in CgCDM 
scenarios differs considerably from the one presented by 
quintesence or ACDM models. As commented in [ p^ , 
the statefinder diagnostic combined with future super- 
novae observations (as, for example, the SNAP mission) 
may be able to discriminate between CgCDM and gen- 
eral quintessence cosmologies. More recently, an analysis 
for the location of the CMB peaks showed that CgCDM 
models and ACDM have very different predictions for 
large values of the parameter As ||l^ . 



ri{z) 



1 



dx 



RoHo J x' f {x T^mi As) 
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where x' = -^^^ = (1 -I- z) ^ is a convenient integration 
variable and the dimensionless function /(x,f2m,^s) is 
given by 



f{x,n,^,As) 



(i-n„,)\ A 



.(8) 



In order to derive the constraints from lensing statis- 
tics in Sec. IV we shall deal with the concept of angular 
diameter distance. For the class of models here inves- 
tigated, the angular diameter distance, Dls{zl, zs) = 
^"11+ zs)^^ , between two objects, for example a lens at 



III. THE AGE OF THE UNIVERSE 

The predicted age of the Universe for the class of 
CgCDM models considered in this paper is given by 



1 



dx 



Ho Jo xf{x,nj^,As 



(10) 



As widely known, a lower bound for this quantity can be 
estimated in a variety of different ways. For instance, Os- 
walt et al. |2ll , analyzing the cooling sequence of white 
dwarf stars found a lower age limit for the galactic disk 
of 9.5 Gyr. Later on, a value of 15.2 ± 3.7 Gyr was also 
determined using radioactive dating of thorium and eu- 
ropium abundances in stars [ p2[ . In this connection, the 
recent age estimate of an extremely metal-poor star in 
the halo of our Galaxy (based on the detection of the 
385.957 nm line of singly ionized ^'^^U) indicated an age 
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FIG. 3: Hoto as a function of the matter density parameter 
for some values of Aa. Horizontal lines correspond to ±2cr 
limits of the age parameter Hoto = 0.95 ± 0.11. 



discussion for the angular diameter distance, for a fixed 
value of Sim the predicted age of the Universe is larger 
for larger values of Ag. If ilm = 0.2 — 0.4, as sugested 
by dynamical estimates on scales up to about 2h~^ Mpc 
II, we find As > 0.96 (see also 



IV. CONSTRAINTS FROM LENSING 
STATISTICS 

In order to constrain the parameters fim and A^ from 
lensing statistics we work with a sample of 867 (z > 1) 
high luminosity optical quasars which includes 5 lensed 
quasars. This sample consists of data from the following 
optical lens surveys: HST Snapshot survey p4|, Cramp- 
ton survey Yee survey |3q, Surdej survey pTt, NOT 
Survey |8) and FKS surveyn|39|. 

The differential probability dr of a beam having a lens- 
ing event in traversing dz^ is ^ 



V RqDos J Ra dzL 



where 



and 



dt 
dzL 



H- 



(i + ZL)/(x,a„,^)' 



dzL, (11) 



(12) 



of 12.5 ± 3 Gyr |2^. Another important way of estimat- 
ing a lower limit to the age of the Universe is dating the 
oldest stars in globular clusters. Such estimates, how- 
ever, have oscillated considerably since the publication 
of the statistical parallax measures done by Hipparcos. 
Initially, some studies implied in a lower limit of 9.5 Gyr 
at 95% confidence level (c.l.) Q. Nevertheless, subse- 
quent studies |2^ using new statistical parallax measures 
and updating some stellar model parameters, found 13.2 
Gyr with a lower limit of 11 Gyr at 95% c.l., as a cor- 
rected mean value for age estimates of globular clusters 
(see also (26|). Such a value implies that the Einstein-de 
Sitter model is ruled out for h > 0.50, while the most 
recent measurements of h point consistently to ft, > 0.65 
p7| , p8| . These results are also in accordance with recent 
estimates based on rather different methods for which the 
ages of the oldest globular clusters in our Galaxy fall on 
the interval 13.8 - 16.3 Gyr |||. 

By assuming to = 13 ± 1 Gyr as a median value for 
the most recent age estimates of globular clusters and 
Ho = 72±8 Km.s~-'^.Mpc~^, in accordance with the final 
results of the Hubble Space Telescope Key Project [p8| , 
we find Hoto — 0.95 ± 0.11, a value that is compatible 
with the estimates discussed above, as well as very close 
to some determinations based on SNe la data |^0[ |3l| . In 
Fig. 3 we show the dimensionless age parameter Hoto as a 
function of Sim for some selected values of As . Horizontal 
dashed lines indicate ±2cr of the age parameter for the 
values of Ho and to considered here. Similarly to the 
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In Eq. (|11|), DoL, Dos and Dls are, respectively, the 
angular diameter distances from the observer to the lens, 
from the observer to the source and between the lens and 
the source. We use the Schechter luminosity function 
with the lens parameters for E/SO galaxies taken from 
Madgwick et al. |||, i.e., = 0.27/1^ IQ-^Mpc^^ a 
-0.5, 7 = 4,?;* = 220km/s and F* = 0.01. 

The total optical depth is obtained by integrating dr 
along the line of sight from zq {z = 0) to zs- One obtains 



r{zs) = ^[Dos{l + ZL)f Rl 



(14) 



In Fig. 4 we show the normalized optical depth as a 
function of the source redshift (zg) for values of Ag = 
0.2, 0.4, 0.6, 0.9 and 1.0. Two different cases are illus- 
trated: a conventional CgCDM model with the matter 
density parameter fixed at Sim = 0.3 and a Chaplygin gas 
+ baryonic matter model with Slf, = 1 — flcg — 0.04. As 
discussed earlier, the reason for considering the later case 
is because one of the strongest claims for a Chaplygin gas 
dark energy is the possibility of a unified explanation for 
the dark matter and dark energy problems ||l^ . In this 
case, one might expect that the only two components of 
the Universe would be the Chaplygin gas and the bary- 
onic matter. Note that in both cases an increase in As 
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FIG. 4: The normalized optical depth (t / F*) as a function 
of the source redshift [zs) for some selected values of As. 
Upper panel: CgCDM models with = 0.3. Lower panel: 
a universe dominated by baryonic matter (p.f, = 0.04) and the 
Chaplygin gas. 



at fixed Vi^ tends to increase tiie optical depth for lens- 
ing. For example, for fim = 0.3, the value of t/F* for 
As = 0.2 at Zs = 3.0 is down from the ACDM {As = 1) 
value by a factor of ~ 2.97, while at the same redshift, 
t/F* for As — 0.4 is down from that for As = 1.0 by ~ 
2.63. By fixing the value of As, for example. As = 0.6, we 
observe that the value of t/F* is smaller for a universe 
with r^ni = 0.3 than for a universe composed only of the 
Chaplygin gas + baryonic matter (17;, = 0.04) by a factor 
of ~ 1.18. This increase of the optical depth as the value 
of As is increased (at a fixed zs and fim) is an expected 
consequence since this model more closely approaching 
the ACDM case as ^ 1. 

The likelihood function is defined by 



C 



Nu 

i=l k=l 



ck^ 



(15) 



where Ni^ is the number of multiple-imaged lensed 
quasars, Njj is the number of unlensed quasars, and p^. 
and p^j. are, respectively, the probability of quasar k to 
be lensed and the configuration probability. These quan- 
tities are defined by 



p^{m,z) 



d{Ae) pc{A0)B{m, z, Mf{A9), Nh) 



B{m, z, Mo, M2) 



(16) 



K B{m,z,Mf{A9),M2) 



PcW 



dzLd{M) 



dzr 



(17) 



(18) 



Af/=Mo(/+l)/(/-l) with /^100-4AmW, (19) 

The magnification bias, B(to, z), is considered in order 
to take into account the increase in the apparent bright- 
ness of a quasar due to lensing which, in turn, increases 
the expected number of lenses in flux limited sample. 
The bias factor for a quasar at redshift z with apparent 
magnitude m is given by p8|, Eg| 



B(m, z) 



B{m, z. Mo, M2) 



(20) 



where 



B{m,z,Mi,M2) 



fdN, 



\ dm 
dNg 



Mo 



Ml 



dM 



(21) 



dm 



(TO + 2.51og(M),z). 



In the above equation [dNq^m, z)/dm) is the measure 
of number of quasars with magnitudes in the interval 
{m,m + dm) at redshift z. Since we are modeling the 
lens by a singular isothermal model profile, Mq = 2, we 
adopt M2 = 10* in the numerical computation. 

For the quasar luminosity function we use Kochanek's 
"best model" [pi 



dNc 



dm 



(to, z) CX (10-'^("-'") + lQ-bi^n-m)yl^ (22) 



where 



Too -|- (z — 1) for z < 1 
m = < Too for 1 < z < 3 

Too — 0.7(z — 3) for z > 3 



(23) 



and we assume a = 1.07 ± 0.07, b = 0.27 ± 0.07 and 
Too = 18.92 ± 0.16 at B magnitude ||. 

Due to selection effects the survey can detect lenses 
with magnification l arg er than a certain magnitude Mf 
given by equation (Q) which becomes the lower limit 
in equation (p2|). To obtain selection function corrected 
probabilities, we follow and divide our sample into 
two parts, namely, the ground based surveys and the HST 
survey. 

FromEq. (|l|) we find that the maximum value of the 
likelihood function is located at fim = 0.4 and As — 1.0. 
At the Icr level, however, almost the entire range of As 
is compatible with the observational data for values of 
ri,n = 0—1. As observed earlier (see Sec. II), this 
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FIG. 5: a) Predicted number of lensed quasars as a function 
of As for Qm = 0.3 (solid line) and Qb ~ 0.04 (dashed-dotted 
line) and image separation A9 < 4. b) Contour for five lensed 
quasars in the parametric space As — ilm- The shadowed 
horizontal region corresponds to the observed range Qrn = 
0.3 ±0.1 H. 

result suggests that a large class of CgCDM scenarios 
is in accordance with the current gravitational lensing 
data. For the sake of comparison, we also analyse some 
possible differences between our best-fit value and the 
one obtained for general quintessence scenarios with an 
equation of state Px = uJxPx 

(XCDM) §. For exam- 
ple, for XCDM models a similar analysis shows that the 
maximum value of the likelihood function is located at 
rim = 0.0 and lo^ = —0.2 Q. Such a model corresponds 
to a decelerated universe with a deceleration parameter 
Qo = 0.2 and a total expanding age of 8.1h~^ Gyr while 
our best-fit CgCDM model corresponds to an accelerat- 
ing scenario with Qo — —0.39 {za = 0.44) and a total age 



of the order of 8.7h^^ Gyr. In Fig. 5a the expected num- 
ber of lensed quasars, riL = J2 Pi (^^^ summation is over 
a given quasar sample), is displayed as a function of Ag. 
As indicated in the figure, the horizontal dashed line indi- 
cates = 5, that is the number of lensed quasars in our 
sample. By this analysis, one finds Ag — 0.9 (ilm — 0.3) 
and As = 0.73 {fit = 0.04). In Fig. 3b we show the 
contour for five lensed quasars in the parametric space 
As — Om. The shadowed horizontal region corresponds 
to the observed range fim = 0.3±0.1 [Q. As a general re- 
sult, this analysis provides flm < 0.45 and Ag > 0.72. We 
also observe that the higher the value of the higher 
the value of As that is required to fit these data. 
V. CONCLUSION 

The search for alternative cosmologies is presently in 
vogue and the leitmotiv is the observational support for 
an accelerated universe provided by the SNe la results. 
In general, such alternative scenarios contain an unkown 
negative-pressure dark component that explains the SNe 
la results and reconciles the inflationary flatness predic- 
tion (r^T = 1) with the dynamical estimates of the quan- 
tity of matter in the Universe (fim — 0.3 ±0.1). In this 
paper we have focused our attention on another dark en- 
ergy candidate: the Chaplygin gas. We showed that the 
predicted age of the Universe in the context of CgCDM 
models is compatible with the most recent age estimates 
of globular clusters for values of — 0.2 and Ag > 0.96. 
We also studied the influence of such a component on the 
statistical properties of gravitational lensing. At la level 
we found that a large class of these scenarios is in agree- 
ment with the current lensing data with the maximum 
of the likelihood function (Eq. ^5|) located at ~ 0.4 
and As = 1.0. As a general result, the predicted number 
of lensed quasars requires Jim < 0.45 and As > 0.72. 
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